Temperature-independent Casimir— Polder forces in arbitrary geometries 
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We show that the Casimir-Polder potential of a particle in an energy eigenstate at nonretarded 
distance from a well-conducting body of arbitrary shape is independent of the environment tem- 
perature. This is true even when the thermal photon numbers at the relevant atomic transition 
energies are large. A compact expression is obtained for the temperature-independent potential, 
which can greatly simplify calculations in nontrivial geometries for experimentally relevant systems 
such as Rydberg atoms and polar molecules. We give criteria for the validity of our temperature- 
independent result. They are illustrated by numerical studies of a particle near a gold sphere or 
inside a gold cylindrical cavity. 
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Introduction. Theoretical and experimental research 
into the Casimir and Casimir-Polder (CP) effects [l] have 
seen a phenomenal surge of activity over the last good 
decade These effects, together with the van der Waals 
forces often collectively termed dispersion forces, are con- 
sequences of ever present fluctuations, quantum and ther- 
mal, of charges, currents and electromagnetic fields. 

Motivated by experimental progress and an increas- 
ing number of technological applications, one has arrived 
at a number of very general insights into the behaviour 
of dispersion forces in recent years. It was shown that 
the Casimir interaction between bodies that are mir- 
ror images of each other is always attractive [3.]- A re- 
cently proven Earnshaw-type theorem generalises this re- 
sult, demonstrating that Casimir forces between dielec- 
tric bodies of arbitrary shapes in vacuum do not support 
stable equilibria General duality symmetry relations 
relate dispersion forces between electric vs magnetic ob- 
jects [5], while scaling relations constrain their depen- 
dence on body separations and sizes Q. 

With experiments typically being conducted at room 
temperature, the impact of thermal photons on the CP 
potential has been of particular practical interest. 
It has recently been studied for various non-equilibrium 
scenarios |8l4ld| where repulsive or spatially oscillating 
(transient) forces have been predicted . On the other 
hand, it has also been shown that the CP potential of 
particles at nonretarded separations (small compared to 
the transition wavelengths) from a metal plate is inde- 
pendent of the ambient temperature, although the avail- 
able thermal photon numbers at the relevant transition 
energies can be very large [13 1. 

Because the geometry and temperature dependencies 
of dispersion forces are closely intertwined 1^, one 



close to conducting bodies of arbitrary shape. The de- 
rived temperature-independent potential is governed by 
the electrostatic Green tensor. For geometries tractable 
via an image-charge method, it is simply the quantum- 
averaged interaction of the atomic dipole with its images 
as produced by the curved body surfaces. 

This very simple solution to an initially complex prob- 
lem of nontrivial geometry out of thermal equilibrium 
could greatly simplify technologically and experimen- 
tally important calculations involving such experimen- 
tally promising systems. It shows that cooling is not 
able to reduce or alter the nonretarded CP potential of 
a particle near any metal body. Two important systems 
for which this temperature-invariance holds in typical ex- 
perimental set-ups are ground state polar molecules 
and atoms in highly excited Rydberg states 
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might suspect this temperature-invariance to be an arte- 
fact of the planar geometry. Quite the contrary, we will 
show in this Letter that this result extends to particles 



these particles, the nonretarded regime extends to sev- 
eral micrometers (polar molecules) and several hundred 
micrometers (Rydberg atoms), respectively. 

We begin by proving that the temperature-dependence 
of the CP force is cancelled at nonretarded separations 
for bodies of arbitrary shapes and give criteria for the 
validity of the temperature-independent result. These 
are illustrated by two numerical calculations for particles 
in nontrivial geometries. 

Temperature-independence in the nonretarded regime. 
In this section we will show the following: Consider a par- 
ticle in the vicinity of metallic bodies of arbitrary shape, 
such that all relevant particle-body separations are much 
smaller than the dominant intra-atomic transition wave- 
lengths (nonretarded regime). Assuming further that the 
reflectivities of the bodies at frequencies up to the atomic 
transition frequencies do not differ from the static reflec- 
tivies, the CP interaction is independent of temperature. 

As shown in Ref. the CP potential of a particle 
prepared in an energy eigenstate \n) may be given as 
Un{r) = U^^ir) + U^{r) with nonresonant and resonant 
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contributions 
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Tr[a„(i^,).ri,.^(r)], 
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which is a good approximation for metals. With these 
two assumptions, the expansion Q for the tensor ^ 
can be approximated by 
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Here, — j2TrkBT/h are the Matsubara frequencies; 
n{(jj) — [e'''^/^'^^^) — 1]^^ is the mean thermal photon 
number of radiation of frequency w at temperature T; 
ctn is the atomic polarisability; and the sum in Eq. ([2]) 
runs over all states \k) to which there exist atomic 
dipole transitions with dipole moment matrix elements 
dnfc = {n\d\k) and energies huJkn = Ek — En- For polar 
molecules 11] and Rydberg atoms alike, only a few 
transitions to neighboring states contribute noticeably to 
the potential. 
The tensor 



r^(r)^^^G«(r,r;c.) (3) 



is given in terms of the scattering part G^^-' of the classi- 
cal Green tensor for the electromagnetic field. We shall 
show that this is a frequency-independent constant in the 
nonretarded limit. To this end it is convenient to use a 
Dyson equation to represent G'^^ as a perturbative ex- 
pansion in xl (1 -I- x/3) where ^) = £(r;a;) — 1 is the 
dielectric susceptibility. As shown in Ref. Il5| . 



G«(r,r';.) = ^^ 

n=0 



x(si;^) 
l + ix(si;w) 



with k ~ ijjjc, pri = Si — r, pi2 = S2 — Si etc. and 



(4) 
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■{[l-ikp~-k^p^]l 



[3 — 3ifc/9 — fc^p^Jepep} 



(5) 



(I: unit matrix; — p/p). We have assumed that all 
bodies are non-magnetic and isotropic, and that r and r' 
lie outside the dielectric bodies, where x = 0- 

In the nonretarded limit \kz\ ^ 1, where z is the 
largest relevant particle-body distance, Eq. ^ reduces 
to 



^(I-3epep)/(47rp3) + 0(fc2) 



(6) 



We further assume that the perturbative parameter is 
frequency-independent in some low-frequency range. 



X(r;w') 



X(r;0) 



l+ix(r;w') i+ix(r;0) 



for \uj'\ < uj, (7) 



X(si;0) 



^ (47r)"+i 

x(s„;0) 



dsi ■ ■ ■ / ds, 



i+ix(si;0) i+ix(s„;0) 



■Ap„. (8) 



for \uj'\ < oj together with Eq. ([5]), meaning T^{r) is 
independent of frequency in this range. 

We apply this result to the CP potentials ([T]) and ([2]) 
in the nonretarded regime |a;fc„z/c| 1. The dominant 
contribution to the Matsubara sum is for < \LUkn\ ^ 
c/z, so we can set Ti^j ~ Tq. Performing the sum accord- 
ing to X)^'o (i^j ) = J2k['^i^kn) + ^](inkdkn/{kBT),we 

obtain 



C/f (r) ~ -1 V[nK„) + i]d„fe-ro(r).d, 



(9) 



Noting that x(r; 0) is real due to the Schwarz refiection 
principle, Eq. ^ shows that Re(r„j^^) ~ Re(ro) = Tq 
in Eq. @. 

The total nonretarded potential hence reads 



-^-y^d„fc-ro(r)-dA; 



(10) 



This remarkably simple expression is manifestly inde- 
pendent of both temperature and transition frequencies 
ujkn- One can explicitly verify that it coincides with the 
zero-temperature result using 27rfcBT'^°^'p — > h d£^ 
in Eq. (P and nriuj) -)> -&{-uJkn) in Eq. ^ [9: unit 
step function]. The temperature- independent CP poten- 
tial is governed by the electrostatic Green tensor Tq. It 
is entirely due to the dipole fluctuations of the particle. 
For geometries where image-charge methods apply, the 
potential coincides with the quantum average of the in- 
teraction energy of the particle's dipole with the image 
dipoles inside the bodies 16|. The potential depends on 
the atomic internal state only via the dipole fluctuations 
X]fc d„fedfe„ = (dd)„. There is hence no qualitative dif- 
ference between forces on ground-state or excited atoms. 
In particular, the sign of the force is the same in both 
cases. 

Note that the temperature-invariance holds for a par- 
ticle in a chosen energy eigenstate or any temperature- 
independent incoherent superposition state. When the 
particle is in a thermal superposition of energy eigen- 
states, the CP potential acquires a weak temperature- 
dependence as discussed in Ref. 13] ■ In addition, the 
internal dynamics of a particle initially prepared in an 
energy eigenstate does of course depend on the ambient 
temperature. 
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Criteria for temperature-independence. The T-inva- 
riant form (ITO|) of the CP potential relies on three condi- 
tions. In the following, we explicitly state these criteria 
and give the leading corrections due to slight violations. 

(A) The relevant particle-body separations must be 
sufficiently nonretarded such that 



ZLOkn 



< 1. 



(11) 



Corrections due to retardation arise from corrections to 
Eq. ([5]). Provided that the approximation ([7]) holds, one 
has r^j ~ To + where the primes indicate w- 

derivatives and x(r; u) — ^ x(r; 0) is understood. We sub- 
stitute this Taylor expansion into Eqs. ([T]) and ([2]). In 
the geometric high-temperature regime kB,T 3> hc/z 3> 
fijwfenl, the nonresonant potential does not contribute, 
because the j = Matsubara term vanishes and all terms 
j > 1 are exponentially small. The retardation correction 
to the CP potential is hence dominated by its resonant 
contribution and it reads 



ctard. 



id„fc-rQ(r)-dA;„. 



(12) 



where rQ(r) is real. It can be reduced by choosing parti- 
cles with larger transition wavelengths or going to smaller 
particle-body separations. 

(B) The reflectivity of the bodies must be sufficiently 
frequency independent such that 



|ro(r)|ri||Rer-':*(r)-ro(r)||«l; (13) 



Here, r"'^°^'(r) is given by Eq. (U) with A approximated 
by Eq. ([6]), but without assuming Eq. (O. Combining 
Eqs. ([9]) and ([2]), the reflectivity correction for fcer ^ 
reads 



d„fe.[RerS-:(r)-ro(r)].dfc„ ^^^^ 



As examples of simple convex, planar and concave ge- 
ometries, consider an atom at distance z from the surface 
of a sphere (radius R <^ z) [17Li_plate [l^l or inside a 
spherical cavity (radius R = z) [18|. In these cases, one 
finds 



AC/, 



retard. 



UniT = 0) 



Cr ctard. 



fceT 

huJhn 



ZUkr. 



(15) 



with Cictard. = -1/3 (sphere), Crctard. = (plate) and 
Crotard. = 3/5 (cavity). The reflectivity corrections read 



rcfl. 



C/«(T = 0) 



fceT z\ujkn\ „ ^ 

-Crefl. T Re 

hWkn C ^/e{uJkn) 



(16) 



with Crcfi. — R/z (sphere), Crefl. = 6 (plate) and Crofl. = 
3 (cavity). For a conductor with plasma frequency ujp 



and relaxation constant 7, the Drude model £{uj) = 1 
1/ 



u!p/[u!{u! + ij)] leads to (w, 7 <C wp, w real) 



Re[i/v/^] ^Up^[^{V^J^+^+\ojm\' . (17) 

The retardation and reflectivity corrections are addi- 
tive to linear order. More generally, the leading correc- 
tion for high temperatures can be obtained by combin- 
ing Eq. ([2]) with the dominant j = contribution from 
Eq. (HI, 

A[/„(r, T)^^Y1 + 0(r°), (18) 



hen 



^kn 



with Ar^(r) = ReFtjj.^ (r) — ro(r); it contains both re- 
tardation and reflectivity corrections. 

As will be demonstrated below, the reflectivity cor- 
rection is geometry-dependent. To wit, Eq. (jl8p distin- 
guishes three classes of geometries where Ar^^ "oc° \uj\^ 
with (i) 77 > 1, (ii) -q — 1 or (iii) 77 < 1. In case (i), the 
reflectivity correction can be made arbitrarily small for 
sufficiently small u, and the T-dependence can be made 
arbitrarily small. For case (iii), the slope of the linear 
correction term oc T reaches a minimum as w decreases; 
it increases when further reducing lo because the reduc- 
tion in the Green tensor difference is overcompensated 
by the growth of the photon number. 

(C) The CP potential must be dominated by transi- 
tions for which (A) and (B) hold. If the conditions hold 
for a transition n — )■ fc, but not for n — > then one must 
have 



«1, 



(19) 



ensuring that the linear T-corrections from retarded but 
suppressed atomic transitions do not contribute. Care 
must be taken in the case where a transition resonates 
with a cavity of quality-factor Q. Then we even require 
that Q(|d„i|/|d„fc|)(fcBr/fi|wi„|) < 1. 

Numerical examples. We now present two illustrative 
examples of non-planar geometries with curvature radii 
on the order of the particle-body separation. For sim- 
plicity, we consider isotropic particles for which 



d„fe-rcj(r)-dfe„ 



j/cptrr^ 



(r). 



(20) 



Consider first the case of a particle at distance r from 
the center of a gold sphere of radius R. The Green tensor 
for this geometry is given in Ref. [20| . In the nonretarded 
limit /cr, kR 1 and for a perfectly conducting sphere 
|£| -> 00, Eqs. (Ho]) and ^ lead to 



C/„(r) 



i?3(6r4 - 3r2i?2 



i?4 



247r£or4(r2 - i?2) 



\dnk 



(21) 



This T-independent potential is in agreement with the 
zero-temperature result found in Ref. |21]. In Fig. [1] we 
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Temperature (K) 



FIG. 1: Comparison of the exact thermal CP potential U{T) 

of a ground-state two-level particle (transition energy hck) FIG. 2: Same as Fig. [T] but for a particle inside a gold cylin- 

outside a gold sphere with the T-independent result Uo from drical cavity. 

Eq. 



demonstrate the accuracy of Eq. (|2ip over a large range 
of temperatures for particles of different transition fre- 
quencies outside a gold sphere. The sphere is described 
by the Drude model with ojp = 9eV and 7 = 35meV. As 
Fig. [1] shows, the correction due to finite temperature is 
about 1% at room temperature when kr = 0.01 and de- 
creases rapidly as fcr — ?• 0, thus exemplifying case (i) with 
•q = 3/2. The temperature correction due to retardation 
(A) is the dominant correction, and the reflectivity cor- 
rection (B) is considerably smaller. 

Consider next a particle inside a cylindrical cavity of 
radius R whose Green tensor can be found in Ref. [lOj. 
In the nonretarded limit and for a perfectly conducting 
sphere, it leads to the potential 

xlliqp)+q'C{qp)}J2\dnk\^ (22) 

(Krm Im- modified Bessel functions). Its reliability is 
illustrated in Fig. [2] In contrast to the previous exam- 
ple, the reflectivity correction (B) dominates; it may be 
shown that except very close to the walls one has case 
(iii) with 77 «i 1/2. Hence as the particle's transition fre- 
quency is decreased, the temperature corrections reach a 
minimum. Note that this problem does not arise when 
suppressing these corrections by reducing the spatial di- 
mensions of the setup. The temperature dependence is 
nonetheless quite modest; for our parameters the linear 
T-correction is at the 5% level at room temperature for 
fci? ~ 0.01 - 0.05. 

Conclusions. We have shown that the CP potential 
of a particle prepared in an energy eigenstate at nonre- 
tarded distance from a conducting body may be approx- 
imated by an expression which is independent of both 
temperature and particle transition frequency, regardless 
of the geometry of the body. In this regime, the interac- 
tion is governed by the transition dipole matrix elements 
of the particle and the electrostatic Green tensor. There 
is no qualitative difference between the potential acting 



on atoms in ground or excited states. For instance, this 
implies that repulsive CP forces that have been predicted 
for ground-state atoms in certain geometries '23] may be 
more easily verified using excited particles in the same 
geometry. 

While the fact that dispersion interactions are of elec- 
trostatic nature in the nonretarded regime is well known, 
the remarkable cancellation of temperature dependent 
contributions only occur once both nonrcsonant and reso- 
nant contributions are fully accounted for. Each of these 
force components, which arise from physically distinct 
phenomena, individually depend strongly on tempera- 
ture and transition frequency. 

The two numerical examples illustrate the reliability 
of the temperature-independent potential. For a particle 
next to a metallic sphere, the temperature dependence 
vanishes for ever smaller values of the retardation param- 
eter ruj/c, whereas for a particle in a cylindrical cavity, 
the T-dependent corrections reach a minimum beyond 
which they cannot be further reduced. The latter effect 
seems to be a peculiarity of the cylindrical geometry. 

This work was supported by the UK Engineering and 
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